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The doped Fullerides can be well described by a Hubbard model, which 
comprises the partly filled, threefold-degenerate t\ u orbital and the on- 
site Coulomb interaction U. The orbital degeneracy is known to shift 
the critical ratio U c /W for the Mott-Hubbard transition towards larger 
values. This puts the half-filled alkali-doped Fullerides A3C60 on the 
metallic side of the transition. Prompted by the recent synthesis of 
isostructural families of integer-doped Fullerides with different fillings, 
we investigate how the orbital degeneracy affects the Mott transition 
at integer fillings n 7^ 3. The calculations are done by fixed- node 
Diffusion Monte-Carlo, using a trial function, which permits us to 
systematically vary the magnetic character of the system. 



1 Introduction 

Due to the relatively large distance between the Ceo molecules in solid Fullerides the 
molecular energy levels merely broaden into narrow subbands. Doping the solid with 
alkali atoms has the effect of filling the subband derived from the LUMO, the t\ u 
orbital. Band structure calculations indicate that the electronic structure changes 
little upon doping, if the compounds retain the same crystal structured The width of 
the t\ u band should thus be fairly independent of doping, and, for the Fullerides with 
fee structure exhibiting orientational disorder, is estimated to be W ~ 0.5 . . . 0.8 eV% 
The Coulomb interaction between two electrons on the same Ceo molecule is in the 
range U ~ 1.2 . . . 1.5 eV, and should also be quite independent of doping! 

An effective Hamiltonian for the electrons in the partially filled band is given 
by a multi-band Hubbard model, taking into account the hopping between the t\ u - 
orbitals on neighboring molecules and the on-site Coulomb interaction U. It turns 
out that, as a consequence of the three-fold degeneracy of the £ lM -orbital, A 3 C 60 is a 
strongly correlated metal, close to a Mott-transition@ Here we investigate the critical 
ratio U c /W at which the metal-insulator transition takes place for the other, integer- 
doped Fullerides A n C 60 with fee structure. Our work has been prompted by the 
recent synthesis of isostructural families of doped FulleridesEI We give an intuitive 
argument for the doping-dependence of the Mott-transition, and present results of 
Quantum Monte-Carlo calculations indicating that the metal-insulator transition is 
shifted towards smaller U c /W as one moves away from half- filling. 



2 Degeneracy Enhancement 



We first consider the effect of orbital degeneracy on the Mott-Hubbard transition in 
a highly idealized model. A system with L electrons becomes insulating when the 
conduction gap 

E g = E(L + 1) - 2 E(L) + E(L - 1) (1) 

opens. In the limit of large Coulomb interaction U ^> t, we can estimate the energies 
entering eqn. (|1|). For integer filling L = nM, where M is the number of molecules in 
the system, all molecules will be occupied by n valence electrons. Hopping is strongly 
suppressed, since it increases the energy by U. Hence for the integer filled system, to 
leading order in t 2 /U, there is no kinetic contribution to the total energy 

E{L) = n(n ~ 1} MU + 0(t 2 /u) (2) 

In contrast, the systems with L = nM ± 1 electrons have an extra electron/hole 
that can hop without an additional cost in Coulomb energy. To estimate the ki- 
netic contribution to the total energy, we calculate the hopping matrix-element for 
hopping of the extra charge to a neighboring molecule. We notice that for the non- 
degenerate Hubbard model a ferromagnetic arrangement of the spins is energetically 
favored (Nagaoka's theorem for bipartite lattices!). For a degenerate Hubbard model, 
however, the hopping matrix element is larger for an antiferromagnetic alignment 
of the spins on neighboring sites. In that illustrated in Fig. [I], instead of 

only one, several equivalent electrons can hop. To estimate the matrix element for 
nearest-neighbor hopping, we calculate the matrix element for hopping of the extra 
electron/hole against a Neel state |0). Denoting the state with the extra charge on 
site % by we find that the second moment of the Hamiltonian (i\H 2 \i) is given by 
the number of hopping channels k times the number of (equivalent) nearest-neighbors 
z times the single-electron hopping matrix-element t squared. By inserting the iden- 
tity in the form J2j where \j) denotes the state with the extra charge hopped 
from site i to the neighboring site j, we then find 

(i\H\j) = Vkt (3) 




Figure 1: Illustration of how an extra electron can hop against an integer-filled background (here: 
degeneracy N = 3, filling n = 3). a) ferromagnetically aligned neighbor: Only one electron can hop. 
b) antiferromagnetically aligned neighbor: There are three hopping channels. 



i.e. the hopping matrix- element is enhanced by a factor of yk compared to the one- 
particle case. Taking the analogy with the one-particle problem further, we then 
expect the energy of the system to be lowered by yk times half the single-electron 
band-width W . The gap in the limit U ^> t is thus given by 



E g = U-(Jk L ^~ 1 + ^T 1 )w/2. (4) 

Without degeneracy we have E g = U — W, while for orbital degeneracy N = 3 we 
find 



filling 


n = 1, 5 


n = 2, 4 


n = 3 


(V^+i + V^-i)/2 


(V2 + l)/2 « 1.21 


(V3 + V2)/2 « 1.57 


\/3 « 1.73 



i.e. extrapolating the result to intermediate U, we expect that the degeneracy enhance- 
ment shifts the critical ratio U c /W to larger values, compared to the non-degenerate 
case. The shift is largest for half-filling and becomes smaller as one moves away 
from n = 3. The symmetry of the result around half-filling reflects the electron-hole 
symmetry (we have implicitly assumed a bipartite lattice). We stress that the above 
argument is not rigorous, its main purpose is to give an idea of how orbital-degeneracy 
affects the electron dynamics. 

3 Quantum Monte-Carlo Results 

To investigate the doping dependence of the Mott-transition in the doped Fullerides, 
we have performed Quantum Monte-Carlo calculations for fee clusters of up to 32 
molecules, using a multi-band Hubbard-Hamiltonian. The hopping integrals have 
been obtained from a tight-binding parameterization! while U has been varied to 
find the Mott-transition. The ground state energies entering in eqn. ([!]) have been de- 
termined by diffusion Monte-Carlo calculations using the fixed-node approximationB 
As trial function we have used the product of a Slater-determinant and a Gutzwiller 
factor. The Slater-determinant has been determined from a Hartree-Fock (HF) cal- 
culation. To allow for different magnetic characters of the trial function, we have 
treated the Coulomb-interaction Uq in the HF calculation as a variational parameter. 

Our results for fillings n — 1 ... 3 are shown in Fig. |2|. As suggested by the de- 
generacy argument, we find that the critical ratio decreases as we move away from 
half-filling. Under the caveat that the neglected effects (higher subbands, multiplet ef- 
fects, electron-phonon coupling) do not change the results drastically, we can conclude 
that the hypothetical doped Fullerides at low temperatures with Fm2>m structure are 
strongly correlated systems with A 2 Cgo being very close to the Mott transition and 
^4iC6o possibly already being in the insulating regime. 




Figure 2: Estimate of the critical ratio U c /W for a multi-band Hubbard model describing the doped 
Fm3m Fullerides at integer fillings n. The bars give our estimates for U c /W. The hatched region 
indicates the U /W range in which these Fullerides are believed to lie. 
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